Scaling Analysis and Renormalisation 

Group for General 
(Quantum) Many Body Systems in the 

Critical Regime 

Manfred Requardt 

Institut fiir Theoretische Physik 
Universitat Gottingen 
Bunsenstrasse 9 
37073 Gottingen Germany 
(E-mail: requardt@theorie.physik.uni-goettingen.de) 



Abstract 

With the help of a smooth scaling and coarse-graining approach of ob- 
servables, developed recently by us in the context of so-called fluctuation 
operators (inspired by prior work of Verbeure et al) we perform a rigor- 
ous renormalisation group analysis of the critical regime. The approach 
is quite general, encompassing classical, quantum, discrete and continuous 
systems, the main thrust going to quantum many body systems. Our cen- 
tral topic is the analysis of the emergent properties of critical systems on 
the intermediate scales and in the scaling limit. To mention some particu- 
larly interesting points, we show that systems typically loose part of their 
quantum character in the scaling limit (vanishing of commutators) and 
we rigorously prove, with the help of the KMS-condition, the emergence 
of the phenomenon of critical slowing down together with the necessity 
of renormalising the time variable. These general features are then il- 
lustrated with the help of an instructive class of models and are related 
to the singular structure of quasi particle excitation modes for vanishing 
energy-momentum. 



1 Introduction 



One of the central ideas of the renormahzation group analysis of, for example, 
the critical regime, is scale invariance of the system in the scaling limit. This is 
the famous scaling hypothesis (as to the underlying working philosophy compare 
any good text book of the subject matter like e.g. PP and references therein). 
Central in this approach is the so-called blockspin transformation, [2]. That is, 
observables are averaged and appropriately renormalized over blocks of increasing 
size. At each intermediate scale a new effective theory is constructed and the art 
consists of choosing (or rather: calculating) the critical scaling exponents, so that 
the sequence of effective theories converge to a (scale invariant) limit theory, 
provided that the start theory lay on the critical submanifold in the (in general 
infinite dimensional) parameter space of theories or Hamiltonians. 

We want to mention a slightly different approach to renormalisation (see for 
example [H] , [01 or [Z| ) , which is more in the spirit of renormalisation in quantum 
field theory. There exist a lot of cross relations but in the following we do not 
discuss these more technical aspects. 

Usually the calculations can only be performed in an approximative way, the 
main tools being of a perturbative character and being typically model depen- 
dent. Frequently, the discussionrelies on spin systems to motivate and explain 
the calculational steps. While the general working philosophy, based on the con- 
cepts of asymptotic scale invariance, correlation length and the like, is the result 
of a deep physical analysis of the phenomena, there is, on the other hand, no 
abundance of both rigorous and model independent results. 

This applies in particular to the control of the convergence of the scaled l- 
point correlation functions to their respective limits if we start from a microscopic 
theory, lying on the critical submanifold. In this case, correlations are typically 
long-ranged and the usual heuristic arguments concerning the manipulations of 
expressions containing block variables of increasing size in the face of long range 
correlations among the blocks becomes rather obscure as one is usually cavalier 
as to the interchange of various limit procedures. One knows from examples, that 
this may be dangerous in such a context. 

Furthermore, the clustering of the higher correlation functions in the various 
channels of phase space may be quite complex and non-uniform. A concise and 
selfcontained discussion of the more general aspects and problems, lurking in the 
background together with a useful series of notes and references, can be found in 
jS], section 7. 

Usually, the crucial scaling relation (the scaling hypothesis) 

WfiLxi, Lxi- /i*) = L-'" ■ L'-^ ■ W^{x^, ...,xi- /i*) (1) 

which is conjectured to hold at the fixed point (denoted by /x* in the parameter 
space), is the starting point (or physical input) of the analysis. Here, Wf denote 
the truncated /-point functions (see below), L is the diameter of the blocks. 
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Lxi denote the centers of the blocks, n is the space dimension, 7 the statistical 
renormalisation exponent. If it is different from n/2 or, rather, the expected 
naive value, we have an ^ anomalous^ scale dimension (for convenience we have 
assumed all observables to scale with the same scale dimension). 

In the following analysis, one of our aims is a rigorous investigation of such 
(and similar) scaling relations for the /-point functions, starting from the under- 
lying microscopic characteristics of the theory. We will do this in a quite general 
manner, that is, the underlying model theory can be classical or quantum, dis- 
crete or continuous. We try to make only very few and transparent assumptions 
as it is our strategy, to deal only with the really characteristic (almost model 
independent) aspects of the subject matter. A central goal of our analysis is a 
rigorous discussion of a number of characteristic properties of both the interme- 
diate and limit states, the observables and dynamics occurring on these levels of 
renormalisation etc., with special emphasis on the quantum aspects. 

In section 2 we develop the conceptual framework and a variety of technical 
tools. As a technical side aspect we discuss the differences between our smooth 
scaling approach and the perhaps more common averaging over sharp blocks. In 
section 3 we show that classical continuum systems can be easily incorporated 
into our general scheme. In section 4 we rigorously study a large class of models 
which can be treated from a unified point of view. We exhibit the close connec- 
tion between the critical exponents and the spectral properties of the correlation 
functions for vanishing energy-momentum. In section 5 we analyse characteris- 
tic properties of the system on the intermediate scales and in the scaling limit. 
Among other things we show that the system may loose some of its quantum 
character in the scaling limit (vanishing of commutators). As a perhaps partic- 
ularly interesting result we provide a rigorous proof of the phenomenon called 
critical slowing down (based on the KMS-condition) together with a renormal- 
isation of the characteristic time scale of the dynamical evolution (see section 
fj5.4|) ). We show, using the class of models studied in section 4, that the dynamic 
scaling exponent, occurring in the renormalisation of the time variable, is closely 
related to the energy- dispersion law of some quasi-particle excitation branch for 
vanishing momentum. 

What regards the general working philosophy, one should perhaps mention 
the framework, expounded in e.g. [4] in the context of the analysis of the ultra- 
violet behavior in algebraic quantum field theory, or, in the classical regime, the 
approach of e.g. Sinai (^19j). While our framework also comprises the classical 
regime (cf. the discussion concerning classical continuous systems in section E}, 
it is mainly designed to deal with the more complicated quantum case. In so far, 
it is an extension of the methods, developed by us in [^, which, on their side, 
were inspired by prior work of Verbeure et al; see the corresponding references 
in jS]. Recently we became aware of a nice treatment of the block spin approach 
in the quantum regime in the bock of Sewell (^SD; who employs methods which 
are different from ours, but are complementing them {quantum (non-) central 
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limits). 

As there exist presumably several thousand papersin this field, we feel unable 
to relate our own approach to all the other approaches or to make a detailed 
analysis of what is entirely new. Our main thrust goes in the direction of a 
conceptual analysis and the developement of a coherent and general point of 
view. In this respect we think, our presentation is reasonably self-contained and 
contains a number of original results. We briefly note that, in order to keep 
the paper within reasonable length, we chose to perform most of the long and 
intricate technical analysis of the scaling behavior of /-point correlation functions 
on the critical submanifold elsewhere. A preliminary treatment of this particular 
problem can be found in the second part of pUj . 

We end this introduction with mentioning a perhaps subtle point. In the 
following we concentrate most of our analysis on the hierarchy of correlation 
functions which can be used to deflne the theory. We generate renormalized limit 
correlation functions from them which happen to be scale invariant (in a sense 
clarified below), thus defining a new limit theory via a reconstruction process. On 
the other hand, we do not openly discuss the flow of, say, the renormalized Hamil- 
tonians through parameter space as a sequence of more and more coarse-grained 
effective Hamiltonians. The characteristics of these renormalised intermediate 
theories are however implicitly given by their hierarchy of correlation functions 
as was already explained in e.g. |H] or 

One should therefore emphasize, that this well-known integrating out or dec- 
imation of degrees of freedom, which characterizes the ordinary approaches is 
automatically contained in our approach! The effective time evolution is carried 
over from the microscopic theory as described in |S| or (in a slightly other con- 
text) in see also P and is redeflned on each intermediate scale, thus implying 
automatically a rescaling or renormalisation of both the time evolution and the 
corresponding Hamiltonian; see section |S| In case we work in an scenario, de- 
flned by ordinary Gibbs states, our framework would exactly yield these effective 
Hamiltonians. Nevertheless, it is an interesting task, to apply our method directly 
to the microscopic Hamiltonian. 

2 The Conceptual Framework 
2.1 Concepts and Tools 

As to the general framework we refer the reader to W^. One of our technical 
tools is a modifled (smoothed) version of averaging (modiflcations of the ordinary 
averaging procedure are also briefly mentioned in the notes in [Hj). Instead of 
averaging over blocks with a sharp cut off, we employ a smoothed averaging with 
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smooth, positive functions of the type 



)1 for T 1 
tor |x| > 1 

Remark: We will see in the following, that the final result is more or less inde- 
pendent of the particular class of averaging functions! 

We note that this class of scaled functions has a much nicer behavior under 
Fourier transformation, as, for example, functions with a sharp cut off, the main 
reason being that the tails are now also scaled. We have 

f^(k) = const ■ i?" ■ f{R ■ k) (3) 

One might perhaps think that this choice of averaging will lead (as a consequence 
of the scaled tails) to a limit theory, being different from one, constructed by 
employing a sharp cut-off. This is however not the case. As the mathematical 
differences between the two approaches, that is, using either sharp or smooth cut 
off functions, are technically a little bit subtle and perhaps not so apparent, we 
discuss some of the technical aspects below. 

We briefly describe the implications coming of translation invariance. We 
have for the correlation functions 

W{xi, .. .,xi) = W{xi - X2, . . . , xi-i - xi) (4) 
The truncated correlation functions are defined inductively as follows (see pi) 
W{x,, xi) = Y,\{ W^{x,,, . . . , (5) 

part Pi 

the sum running over all partitions of the set {xi, . . . The (distributional) 
Fourier transform reads 

iy^(pi, . . . ,J9;) = W'^{pi,Pl+P2, • • • ,Pi H Pi-i) ■ 8{Pi H Pi) (6) 

The dual sets of variables are 

i 

yi := Xi - Xi+i , Qi = "^Pj «<(/-!) (7) 

i=i 

In contrast to the averaging procedure, introduced above, the usual block- 
variable-SLveiSigmg is a sharp cutoff averaging, performed for example over balls, 
Bpi, of radius R. That is, observables are integrated over balls, Bji, with the help 
of the incidence functions 

fl for \x\ *^ 1 
n ^ I 1 
tor \x\ > 1 
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The averaging over operators, leading to the so-called block or fluctuation oper- 
ators, is performed in the following way: 



with the exponent 7 suitably chosen and const being a possible numerical and 
unimportant multiplicative constant (related e.g. to the volume of the unit ball 
or something like that). A corresponding expression holds for xr replacing f^. 

Remark: Here and in the following, the A''s are always normahzed to {A) — 0, in 
order to really get the pure fluctuation effects. 

At this point we want to state a general principle which allows to choose an 
appropriate scaling exponent, 7. As in most of the discussions in the literature 
only a particular fixed field, (f){x), or spin, S{xi), is employed, it is frequently not 
clear that something has actually to be said in this context. This holds the more 
so as quite a few different renormalisation schemes (or philosophies) are used in 
practice, with the tacit understanding that the critical exponents are physically 
apriori given and insensitive to the concrete decimation procedure. This problem 
becomes, in our view, more virulent in the quantum regime with, usually, a lot 
of different non-commuting observables. 

We think that, if we adopt a more general viewpoint, the necessary general 
principles become more apparent. This holds also for what we call a possible 
problem of consistency. This problem consists of the following points. 

Observation 2.1 

1. It is reasonable to choose the scaling exponents so that certain two-point 
auto-correlation functions of block observables survive in the scaling limit. 
Note that an inappropriate choice drives the auto-correlation functions ei- 
ther to zero or infinity! 

2. The Cauchy-Schwarz inequality then guarantees that at least the corre- 
sponding mixed two-point functions remain finite in the limit. 

3. On the other hand, this shows that one may have a certain freedom in 
selecting the correlation functions and, by the same token, the observables 
one wants to survive in the limit. This will of course affect the structure of 
the possible limit theories. 

4. If, on the other hand, we have a lot of different (non-commuting) quan- 
tum observables together with their composites, it is presumably not an 
easy task to make all these (possibly independent) choices in a consistent 
way so that a coherent Hilbert space structure results in the limit. This 
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problem becomes virulent if we end up with a theory having non- vanishing 
higher truncated correlation functions. The reason is that possible obstruc- 
tions may result from the decay behavior of higher n-point functions in the 
difference variables which has to be in complete balance with the chosen 
scaling exponents. These, on their side, are already fixed by the 2-point 
functions! We briefly discussed this issue in the last section of [H] and we 
make a more detailed analysis in the second part of our paper jHUj . 

Conclusion 2.2 We fix the renormalisation exponents, 7^, of the respective 
observables via the non-vanishing and finiteness of (a class of) 2-point auto- 
correlation functions. This will yield constraints on the scaling behavior of higher 
correlation functions, the consistency of which we then can check. 

In the following we will mainly employ the smooth cut-off procedure which 
leads to a more transparent behavior of various expressions in Fourier-space. It 
is satisfying that in the cases, we can actually control, it leads to results being 
identical to the version with sharp volume cut-offs. In order to compare these 
two cut-off conventions we study in a first step various peculiar properties of 
the averaging functions, Xr{^)- The Fourier transform of the smooth functions, 
fnix), are again smooth, living in the Schwartz-space, S, i.e., decrease fast to- 
gether with all their derivatives. In jH] we crucially employed or properties 
of various expressions. In contrast to fR{k), the Xi?(^)'s are no longer in as 
Xr{x) has a jump discontinuity. On the other hand, it is in as 

00 > I \xix)R\'rx = I \xR{k)\'d''k (10) 

We have the little lemma 

Lemma 2.3 The Fourier transform, XR{k), is in C°° fl Cq but not in L^. It is 
however in L^. We have the same scaling behavior for X-r(^) as for fR{k), that 
is 

XR{k) = const ■ ■ Xi{R ■ k) (11) 

Proof: The first statement follows from the compact support of Xr{x) and the 
Riemann-Lebesgue lemma. The second statement follows as in the smooth case. 
□ 

An explicit calculation for n = 3 yields: 

Ak\ 

Xi{k) = const ■ \k\^^ ■ I r-sm{r)dr (12) 

Jo 

For |A;| ^ the integral is proportional to \k\^. Furthermore we can show that 
the expression is in fact infinitely differentiable in \k\ =0. For ^ 00 a partial 
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integration yields an expression proportional to — \k\ ■ cos{\k\) +sin(|A;|). That is, 
we have in leading order for \k\ oo: 

Xi(A;) ~ |A;|-2 forn = 3 (13) 

We mention some further peculiar properties of the indicator function, Xb{x), 
not shown by other functions. From Xb{x) = Xb{x) ■ Xb{x) we infer for the 
Fourier transform 

XB{k) = Xb * XBik) (14) 
and correspondingly for higher powers. 

Corollary 2.4 By Young's inequality (see e.g. jSHj) we know, that in general 
the convolution of L^-functions is only in L°°. The preceding formula shows that 
the convolution of Xsik) with itself is again in L^. 

Note that such a result is not immediately evident from the concrete form of 
the respective Fourier transforms. In the case n = 1 say, the Fourier transform 
is essentially of the form sin(A;)/A;. The result for the convolution comes about 
due to the peculiar oscillatory character of the expression and would not hold for 
e.g. |xi(/c)|. We will briefly analyse in the following subsection to what extent 
the renormalisation process is influenced by these slightly nasty features of sharp 
cut-off functions. 



2.2 The case of Normal Fluctuations 

As in 1^, we assume that away from the critical point the truncated /-point 
functions are integrable, i.e. G L^{R^^^~^^)^ in the difference variables. 
Hi := Xi — Xj+i. As observables we choose the translates 

AR{ai), Aniai) , A^(a) := i?-"/^ . j Aix + a)f{x/R)rx (15) 

(where, for convenience, the labels 1 . . .1 denote also possibly different observ- 
ables). We then get (for the calculational details see [S], the hat denotes Fourier 
transform, translation invariance is assumed throughout, the const may change 
during the calculation but contains only uninteresting numerical factors): 

{An{a^)---AR{ai)f = const ■R'''/^- 
J fiRpi) ■ ■■fi-R[pi + - ■ ■+pi-i])-W^ipu . . . ,pi.,)-e-'^'"P"'^-e''''^""P'lldpi 

= const ■ . ^-(i-l)n. 

/ /(K) • ■ ■ fi-[p[+- ■ ■+pU]yw^{p[/R, . . . ,p;_i//?)-e-^s"^(^^/^)--e-' JJ 
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We now scale the a^'s like 



ai:=R-Xi, Xi fixed (17) 

This yields 

{AR{R-X,)---Aii{R-Xi)f = 

const ■ i?(2-/)n/2 . [ g-.E'rV:x, . ^^x, E'r^L 



fip'i) ■ ■ ■ /(-[K + ■ ■ ■ + P'l-i]) ■ w^ip'i/R, ■ ■ -M-i/R) n dp'^ (18) 

As the / are of strong decrease and continuous and bounded by assump- 
tion {W^ G L^(R"('~^))!), we can perform the limit i? — oo under the integral 
(Lebesgues' theorem of dominated convergence) and get: 



Case 1 (/ > 3): 



Case 2 (/ = 2): 



hm {Ar{R ■ Xi) ■ ■ ■ Ar{R ■ Xi)f = (19) 

it— >oo 



lim {Ar{R ■ X,)Ar{R ■ X2)f = const ■ [ W^^(O) ■ e-^^'i^^^"^^) ■ f{p[) ■ H'P'M 

(20) 

We arrive at the conclusion 

Conclusion 2.5 Assuming L^-clustering in the normal regime away from the 
critical point and employing a smooth cut-off, all the truncated correlation func- 
tions vanish in the limit R ^ oo apart from the 2-point function. We hence have 
a quasi free theory in the limit as described in |0 or in the work of Verbeure et 
al (cf. the references in P) 

In the case of smooth averaging we employ the transparent behavior of the 
Fourier transformed expressions. On the other hand, the Fourier transform is in- 
herently non-local, which sometimes makes the analysis more complicated. When 
using instead the sharp cut-off convention, we described above, the behavior of 
the respective Fourier transforms becomes opaque in the general case. On the 
other hand, we can try to stay in coordinate space and perform the analysis 
there. Proceeding as in the smooth case but avoiding Fourier transformation we 
get after some straightforward manipulations 

(AniRX,) . . . An{RXi)f = ^"'"Z' " " / W^y'i, • • • , vli) 
l-l l-l 
■ xiiJ2 ^''y'j + ^3 -Xi)--- xi{x'l - Xi)dy[ . . . dx'l (21) 

i=i i=i 
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with Xi — xi — ^ Uj by the definition in the preceding subsection and 

i-i 

x',^Xi + RXi, x',-x'i^J2y'j^ <--AlR (22) 

i 

Again, the hmit can be performed under the integral and is zero for I > 2. 
For the two-point function we get 

^mJAniRX^) • AR{RX^)f = j W^{y)dy ■ J Xi{x - Y)xi{x)dx (23) 

with Y := Xi — X2. We hence conclude: 

Conclusion 2.6 In the normal situation of L^-clustering of correlation functions 
and sharp cut-off functions, Xr{x) = Xii^/R), we get the same results as in the 
case of smooth cut-off functions. However, to prove this, we have to perform the 
analysis in real space and avoid Fourier transformation. 

Corollary 2.7 It is obvious from the preceding discussion that the particular 
form of the averaging functions, /r (x) = f{x/R), need not even simulate a volume 
averaging. For the argument to hold, it is e.g. sufficient that f{x) is bounded 
with /(O) 7^ and has compact support. What only changes is an unimportant 
multiplicative factor, depending on the type of function, being chosen. 

2.3 The Relation to the Heuristic Scahng 
Hypothesis 

In the following sections we develop a rigorous approach to block -spin renormal- 
isation in the realm of quantum statistical mechanics, which tries to implement 
the physically well-motivated but, nevertheless, to some extent heuristic scaling 
hypothesis. The analysis will be performed both in coordinate space and Fourier 
space. In this subsection we restrict our discussion to the two-point correlation 
function, for which the asymptotic behavior is simpler and more transparent. 

Remark: In the rest of the paper we replace the exponent n/2 in the definition 
of AR{a) by a scaling exponent 7', which will usually be fixed during or at the 
end of a calculation. It plays the role of a critical scaling exponent and is model 
dependent. 

Let us hence study the behavior of 

{Ar{R ■ X^)Ar{R ■ X2)f = R-^^' ■ J W'^iixi - X2) + R{X^ - X2)) 

• f{xi/R)f{x2/R)dxidx2 

= i?-27'+2n J W^'iRiix, - X2) + (Xi - X2)]) ■ f{xi)f{x2)dXidX2 (24) 
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We make the physically well motivated assumption that, in the critical regime, 
W'^ decays asymptotically like some inverse power, i.e. 

W^{xi-X2) ~ {const+F{xi-X2))-\xi-X2\'^'"~''^ 0<a<n , F{x) G (25) 

for \xi — X2I —>■ 00, F bounded and well-behaved. 

From the last line of we see that, as / has compact support, we can 
replace W^, for (Xi — X2) 7^ and i? ^ 00 by its asymptotic expression and get 
for R large: 

{Ar{R ■ X,)Ar{R ■ X2)f ^ const ■ i?-27'+2n . ^-(n-a) . J |^ _^ y | -(n-a) . f ^ j ^y^^y 

(26) 

We choose now 

7' = (n + a)/2 (27) 
and get a limiting behavior (for R —>■ 00) as 

const- J \y + Y\~^-"^^-f*fiy)dy (28) 

with y = xi — X2, Y = Yi — Y2 and 

f*f{y) ■■= I f{y + X2)- f{x2)dX2 (29) 



We see that in contrast to some of the general folklore, the limit correlation 
functions are not automatically strictly scale invariant but still depend in the 
above integrated manner on the chosen smearing functions, /. Full scale invari- 
ance is recovered in the regime Y 00. Central in the renormalisation group 
idea is that systems on the critical surface (i.e., critical systems) are driven to- 
wards a fixed point, representing a scale invariant theory. This idea is usually 
formulated in an abstract parameter space of, for instance, Hamiltonians. In our 
correlation function approach scale invariance at the presumed fixed point would 
prove its existence via the scahng properties of the correlation functions, that is 

W^{L • (X - F); fi*) = L-2(n-7')|yT(^ _ y. (30) 

with fi* describing the fixed point in the (usually) infinite dimensional parameter 
space. We see from the above that this picture is asymptotically implemented 
by our above limiting correlation functions, as we have (with the choice 7 = 
(n + a)/2): 

iyj(x-r;^*) ~ (31) 

in the asymptotic regime. That is, the above scaling limit leads to a limit (i.e. 
fixed point) theory, reproducing the asymptotic behavior of the original (micro- 
scopic) theory. 

One should however note that in the more general situation of /-point cor- 
relation functions we have to expect a more complex decay behavior and the 
existence of various channels as varying clusters of observables move to infinity. 
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2.4 Strategies for a Renormalisation Analysis on the Crit- 
ical Surface 

Typically, the numerical scaling analysis is developed for the system being away 
from the critical surface. The reason is that away from criticality, under the 
heuristic assumption of e.g. exponential clustering, the analysis is not beset 
with technical difficulties as, for example, the interchange of limits and dealing 
correctly with long range tales in correlation functions. It is then frequently 
argued that, in case the system is sufficiently near to the critical surface, the 
orbits of renormalized model systems nevertheless will approach the vicinity of 
the fixed point, so that one can make a linear stability analysis of eigenvalues of 
the renormalisation group around the fixed point. The philosophy is that these 
systems will ultimately leave the vicinity of the fixed point. 

In the second part of IHOj (see also the last section of 8J) we undertook to 
sketch a rigorous renormalisation framework for systems, lying on! the critical 
surface. Due to the inherent long-range correlations, one must be extremely care- 
ful in performing such an analysis. As such a rigorous analysis is both technically 
demanding and a little bit tedious and incorporates a variety of interesting math- 
ematical side aspects like e.g. a singularity analysis of distributions and pseudo 
differential operators, we decided to separate this rather technical investigation 
off and give only a brief discussion of one of the methods in this subsection, which 
we exemplify with the help of the pair correlation function. 

The general idea is it, to extract and isolate the characteristic singular be- 
havior of the correlation functions which is responsible for the weak decay of 
correlation. With W'^{x) the truncated two-point function, we, making the pre- 
ceding analysis more rigorous, assume the existence of a certain exponent, a, so 
that (x^ denoting the vector-norm squared) we can make the following decompo- 
sition. 

G{x) := W^ix) ■ (1 + x2)("-")/2 _ ^^^^^ ^ (32) 

with a decaying (non-singular) F which is assumed to be in L^. Fourier transfor- 
mation then yields: 

R-^^ ■ j W^{{xi - X2) + R{Xi - X2))f{xi/R)f{x2/R)dxidx2 

= R-^^ ■ j G((xi - X2) + i?(Xi - X2)) ■ [1 + ((xi - X2) + R{X, - X2))2]-("-")/2. 

/ (xi /-R) / (x2 / R) dxidx2 

J ^_,H,(,,_,,)^^_2 ^ ((^^ _ ^ _ X2))2)-("-°)/V(a;i)/(x2)rfXidX2 

(33) 
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where we made the substitution x ^ R - x. 

As the support of / is in principle arbitrary, we now assume it to be contained 
in a sufficiently small ball around zero (or, alternatively, {Xi — X2) sufficiently 
large so that (xi — X2) + {Xi — X2) 7^ for Xi in the support of /). With 

G{p) = const ■ 5{p) + F{p) (34) 

the leading part in the scaling limit — 00 is the 5-term. Asymptotically we 
hence get for _R ^ 00 (setting y := xi — X2 Y := Xi — X2) the result, already 
conjectured in the preceding subsection: 

. ^^^^^ ■ j \y + ^^r^"""^ ■ / * fiy)dy (35) 

with 

f*f{y)--= I f{y + X2)-f{x2)dX2 (36) 



and y + Y on supp{f). 

The reason why the contribution, coming from F{p), can be neglected for R 
00 is the following: / is assumed to be in P; by assumption the prefactor never 
vanishes on the support of f{xi). Hence the whole integrand in the expression in 
square brackets is again in V and therefore its Fourier transform, g{p'), is in S 
(with p' := Rp), that is, of rapid decrease. We can therefore perform the i?-limit 
under the integral and get a rapid vanishing of the corresponding contribution in 
i? for i? — > 00. 

lim R-n ■ I F{p'/R) ■ e"*^'^ ■ gip)dy = (37) 

R^co J 

As / * / has again a compact support, we have that, choosing 

7 = (n + a)/2 (38) 

the limit correlation function behaves as ~ X2|~^"~"'' as in the more heuristic 
analysis heuristic analysis of the preceding subsection. 

Along these lines, or choosing a slightly different method (see ISO]); one can 
proceed in the more difficult case of /-point functions. In the course of this 
analysis an interesting phenomenon does pop up which leads to some remarkable 
constraints as to the consistency of the whole renormalisation picture. The critical 
exponents are typically fixed by the assumed non-vanishing of the scaled auto- 
correlation functions. On the other hand, the truncated /-point functions may 
have a much more intricate cluster behavior (having, in particular, a variety 
of decay channels). If one wants to go beyond quasi- free limit theories, some 
higher truncated correlation functions have to be non-vanishing in the scaling 
limit. For this to be the case, there has to be some fine-tuning between their 
decay behavior and the values of the critical exponents, which one got from the 
two-point functions. 
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3 Some Remarks on Classical Statistical Sys- 
tems 



In this section we want to briefly indicate how our framework can be implemented 
in the regime of classical statistical mechanics. The situation is more or less 
obvious in the class of spin- or lattice-systems. The translation group is replaced 
by some discrete lattice group. The Fourier vectors run through some Brillouin 
zone instead of M" , while in coordinate space we employ the same kind of smearing 
and averaging functions as in the continuous case, the only difference being the 
replacement of integrals by sums. For continuous classical KMS-systems, some 
more words are perhaps in order (cf. e.g. [20], |22|; 

As infinitely extended phase space , X, we take the set of sequences, x, 

x={r,^Pi)T=i = {xi)T=i (39) 

of points 

x^ = (r„ Pi) e X M" , r, ^ r, for i ^ j (40) 

having the local finiteness property, i.e., the number of points, r^, occurring in x, 
is finite in each bounded set , V G M"". 

As local m-particle observables, A^"^\ we take 

A(-\x):= f^"'\xn.....x,J (41) 

ii<---<i77i 

/ from the class of smooth function with compact support in coordinate space 
(the details are of course a matter of convenience). Poisson brackets can then be 
defined as usual and the local finiteness property guarantees that the expression 

oo 

{A, B}{x) := J2i9A/drj ■ dB/dpj - dA/dpj ■ dB/drj) (42) 
j=i 

is well-defined. 

The thermodynamic equilibrium states are now probability measures on the 
Borel-(T-field, defined on the phase space equipped with the topology canonically 
induced by the class of observables. The classical KMS- condition we usually 
employ in the form [A, B real): 

{{A,B})=^{B{A,H}) (43) 

H being the Hamiltonian. 

As in quantum statistical mechanics, we can define certain distributional point 
fields or densities at, say, coordinate r over the phase space, like e.g. 
particle density: 

nr{x) := ^ 6{r - Vi) (44) 
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momentum density: 

Pr{x) ■.= '^Pi- S{r -Ti) (45) 

i 

energy density, stress tensor density etc. (see in particular and [25 where 
these notions have been systematically employed). 

The l-point distribution functions can hence be expressed as follows: 

p«(ri,...,n) := 5^ (5(ri-riJ---5(n-rO) (46) 

il<-<il 

iu running through the indices occurring in x. Ordinary observables can be 
reconstructed by integrating these densities over local test functions. For a one- 
particle observable we have for example: 

Af := j a{r) ■ f{r)d''r (47) 

with a(r) a one-particle density and correspondingly for more complex densities. 

From these remarks one sees immediately, that the whole procedure, we de- 
velop in the following, can be immediately transferred to the regime of classical 
statistical mechanice without significant changes. 



4 A Class of Examples 

We argued that in the case of poor, that is, non-integrable clustering, it appears 
to be mathematically more reasonable to perform most of the necessary analysis 
in coordinate space, as the behavior in Fourier space my be quite involved in the 
vicinity of (cu, k) = (0, 0). 

The situation improves however if one has a more precise knowledge of the 
form of correlation functions in Fourier space near (0, 0). We note in passing that 
our approach is by no means restricted to the case of critical systems. It does also 
apply to systems at zero temperature or systems above or below a phase tran- 
sition line. One may have more precise information in Fourier space in various 
situations like e.g. spontaneous symmetry breaking (see fl^ and [H] and further 
literature given there) or for particular correlation functions and/or commutators 
(so-called sum-rules) . The relevant contribution in e.g. the 2-point function can 
stem from sharp excitation branches or excitations, having a finite lifetime, which 
is the typical situation in interacting many-body systems. 

Remark: The assumptions in the following discussion can be considerably weak- 
ened and are only made to cover a sufficiently general and coherent class of 
models. 

In order to better understand the effects of our general scaling approach, we 
deal in this section with a fairly large class of relatively manageable and simple 
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models at non-zero temperature which belong to the group of quasifree systems. 
Note however that in contrast to, say, relativistic quantum field theory, we have 
in general no strong covariance properties. That is, even quasifree systems are 
not completely uninteresting and supply us with a whole bunch of useful model 
systems approximating important non-trivially interacting systems. As this no- 
tion slightly varies from author to author, we make the following assumptions. 

Assumption 4.1 Our class of models is assumed to have the following properties 
(in addition to the usual standing assumptions, we do not repeat here; see e.g. 

m) 

The KMS-representations, np, of the quasi-local algebra {A, at^x) are assumed to 
be quasifree and faithful, that is 

1. All n-point functions are products of 2-point functions. 

2. TV/siA) 7^ if A ^ in ^. 

The second assumption seems to be physically reasonable (and can in principle 
be weakened) as it avoids redundancies but need not! be fulfilled in general. 
This situation occurs of course when the original algebra has a non-trivial center 
and one studies representations which are factors, in which central elements are 
mapped onto c-numbers. Note that tu/j, the KMS-state, is always faithful (that 
is, separating) in 7r^(^)" (the GNS-representation) , however this need not be the 
case with respect to A itself. 

This point is relatively subtle from a more physical point of view and not much 
seems to be known. There is a discussion in jTn],p.85ff. which is based on the weak 
closure. A" of the original algebra. But, typically, an equilibrium state is given 
via its local restrictions in form of Gibbs-states, that is, it is naturally only defined 
on quasi-local elements of the algebra and not on the weak closure. If A is simple, 
the representation is faithful. (For an example of a non-faithful representation 
see |2ni)- Note that the so-called order parameters, the non- vanishing of which 
usually signal the occurrence of new phases, are typically global "observables" 
(for example, meanvalues, not belonging to the algebra of quasi-local observables) 
and are c-numbers in pure phases, i.e. factor states. 

The above assumptions have both a simple technical consequence and a con- 
sequence which is perhaps remarkable from a more physical point of view. 

Lemma 4.2 Under the assumptions being made the commutators in each rep- 
resentation, TTp, are c-numbers which do not depend on the KMS-state, that is, 
in contrast to the 2-point functions, they are state-independent. 

Proof: i) The c-number property follows immediately from the vanishing of all 
higher truncated correlation functions and is in fact independent of the other 
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assumptions. With span (77/3(^)^/3) being dense in the GNS-Hilbert space and 

(0/3,71/3(71) • [7r/3(5(a;,t)),7r/3(C)] -77/3(^0^^/3) (48) 
being a sum of 2-point functions, this expression can be shown to be equal to 

(f]/3, [7Tf,{B{x,t)),7lf,{C)]nf,) ■ (l^/3,7r/3(A) • 77/3 (A') 1^/3) (49) 

ii) The faithfulness of implies (with CQQ{x,t) a function, which follows from 
i)) 

C^Bci^,t) = MB{x,t)),Trp{C)]=7r^{[B{x,t),C]) = CBcix,t) = [B{x,t),C] 

(50) 

with cbc{x: t) — C!qq{x, t) being independent of the concrete KMS-representation 
as 

7r/3([i?(a;,t),C]-C^c(^,^)-l) = (51) 

□ 

The physical relevance of the above observation is the following. With 

F^s(x,t) := (Q^,7r/3(A)(x,0 -71^(5)0/3)^ (52) 

and 

Cij,{x,t) := (Q^, [Kp{A){x,t),7:p{B)]np) (53) 
we have the general expression for the respective Fourier transforms 

F^^(a;, k)^{l- exp(-/3a;))-i • Ci^{uj, k) (54) 

Usually, both F and C depend on the parameters, fixing the KMS-state. Our 
assumptions guarantee that for our model class the temperature dependence on 
the rhs is entirely concentrated in the prefactor, (1 — exp(— /9a;))"^, that is, we 
have 

CorollEiry 4.3 For our model class it holds 

F^5(a;, k)^{l- exp(-/g^))-i • cab{uj, k) (55) 
with CAsi'^-ik) temperature independent. 

As cab{^-, k) is universal, it is typically simple to calculate; use e.g. some ground 
state representation. 

For the further analysis we choose A,B selfadjoint and get for C{u!,k) (we 
supress the labels A,B): 

C{uj, k) = F{uj, k) - F{-uj, -k) = (1 - exp(-/5a;)) • F(a;, k) (56) 
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and hence 

Re F{-uj, -k) = exp{-(3iu) ■ Re F{uj, k) (57) 



Im F{-uj, -k) = - exp{-f3uj) ■ Im F{uj, k) (58) 

thus clearly exhibiting the two-sidedness of the (tu, /c)-spectrum in temperature 
states. 

As, in contrast to the relativistic context (cf. e.g. [22] or jSHj), we have in 
general no strong covariance and/or spectrum conditions for the 2-point functions, 
we have to make some reasonable assumptions which are fulfilled in typical many- 
body systems (for more details see [23] and |2H])- 

Assumption 4.4 We assume that the excitation spectrum of -^(0;, k) fulfills 
F{u!,k) = F{u!,—k) and contains a sharp excitation branch {e{k) = e(|fc|)), 
describing stable quasi particles or collective excitations, with the remaining part 
being integrable and absolutely continuous around (cu, k) = (0, 0). We denote the 
singular contribution by 

F^ingiuJ, k) ■= 4{k) ■ 5{iu - {e{k) - /i)) + J^{k) ■ 5{iu + ie{k) - fi)) (59) 

Remark: Note that in the translation invarinant case the above Fourier transforms 
are measures! 

From the above relations we conclude that 

Re J^_{k) = exp (-/5(e(fc) - jj,)) ■ Re J^{k) (60) 

Im J^(fc) = - exp {-f3{e{k) - jj,)) ■ Im J^{k) (61) 

with fi the (temperature dependent; in case temperature and density are chosen 
as independent parameters) chemical potential. We arrive at the following result: 

Lemma 4.5 We have 

J^{k) = exp (-/3(e(A;) - /i)) ■ J^{k) (62) 

and 

C^ingioo, A;) = (1 - exp (-/3(e(A;) - /i))) • 4{k)6{uj - (e(A;) - /x)) 

- (1 - exp (-/3(e(A;) - /i))) • 4{k)6iuj + (e(A;) - /i)) (63) 

As Csingi^, k) has to be independent of (3 for our class of models, we have fur- 
thermore 

Jl{k) = (1 - exp (-/3(e(A;) - ^)))-i ■ j^k) (64) 
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Proof: This follows directly from the preceding formulas. □ 

As our commutator function is universal, it should not contain the typical 
singularities which show up in connection with phase transitions and critical 
phenomena. As to this point we refer to the discussion in e.g. and |28j . 
These phenomena are typically representation dependent. Therefore, on physi- 
cal grounds, the function j{k) should be bounded near k = and e{k) can be 
identified with the dispersion law of an elementary excitation which, in the non- 
relativistic context, for short-range interactions, passes through zero for k —>■ 0. 

Remark: In |25j we discussed various dispersion laws. Frequently a simple power 
law behavior prevails. 

We now apply our scaling procedure to the class of model systems described 
above. In a first step we want to choose the scaling exponent, 7 = 7^, in the 
expression 

Ar = R-^ ■ j A{x + RX) ■ fR{x)(rx (65) 

so that the corresponding autocorrelation function (remember the standing as- 
sumption {A) = 0) 

{Ar{RX,) ■ AniRX^)) (66) 

is both finite and non- vanishing in the limit R 00. 

Our above made observations or assumptions about the spectrum of the 2- 
point functions show that, provided we have a more detailed knowledge of the 
system under discussion, we can, even in the case of long-range correlations, 
perform the analysis in Fourier space getting 

{Ar{RXi) ■ Ar{RX2)) = 

const ■ i?-27+" . j exp(-zA:'(Xi - X2)) ■ f{k')f{-k') ■ Faa{^, k' / R)dujdk' (67) 

with k' := Rki and f{k')f{—k') = |/(A;')p for f{x) symmetric and real. 

In the following we are concerned with the renormalisation of the singular part 
of the spectral contribution as the absolutely continuous part is (by assumption) 
harmless. For u > we have to consider the term 

lim Ir := 

lim i?-2(7-n/2). f exp{ik{X^-X2))il-exp{-P{e{k/ R)-fx)))-^-j{k/ R)-\f{k)\^d''k 

-R— >oo J 

(68) 

We do not intend to discuss the mathematically most general case but rather 
concentrate on situations which are reasonable from a physical 
point of view. 
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Assumption 4.6 We assume that in leading order e{k) behaves hke 

e{k) - |A;|" for ^ (69) 

with a > 0. 

There is the possibihty that j{0) is finite but non- vanishing or that j{k) vanishes 
for /c — > 0. We begin with the discussion of the case of non- vanishing j{0). 

I) e{k) ~ |A;|" near k = 0, j{k) continuous and 7^ in A; = 

Observation 4.7 i) For 7^ nothing pecuhar happens and we are in the 
normal situation with 7 = n/2. 

ii) For /X = 0, the typical situation at or below the critical point, we have the 
following behavior 

Iji^ R-^i^-i^+°^y^) for R^oo (70) 

hence, the anomaleous scaling dimension is 

7=(n + a)/2 (71) 

with a the exponent in the dispersion law of the sharp elementary excitation 
mode. 

II) j{k) vanishing in A; = 

This situation is by no means entirely exceptional. Take for example the time 
derivative at t = of the observable A. In the spectrum of the autocorrelation 
function this leads to an additional prefactor, oj^, in front of FAA{uj^k). In the 
singular contribution, J+{k), this results in an additional factor, e(A;)^, and hence 
in an additional contribution in the scaling exponent 

79A = J^/S + q;/2 - a = 7^ - a (72) 



Observation 4.8 If one wants dtA to be a non- vanishing observable in the scal- 
ing limit, its scale dimension has to be chosen as 

IdtA ^-fA-a (73) 

Similar considerations have to be made for other functions of elementary ob- 
servables. If the spectrum is known qualitatively as in our case, this can in fact 
be done in every concrete case. Note furthermore that the temperature indepen- 
dence of the commutator is technically convenient but not absolutely necessary. 
The same conclusions do hold if the spectral weight along the sharp excitation 
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branch is temperature dependent. However, in that case we do not have an apri- 
ori knowledge as to its precise form which may vary with (3. One can also treat 
the case of excitations having a finite lifetime (cf. jlH])- The excitation branch 
now has a finite width and the calculations become even more model dependent. 
On the other side we proved in [2H] that for /? 7^ sharp excitation branches 
typically belong to elementary excitations having no interaction with the rest of 
the system. 

Remark 4.9 As to interesting consequences concerning the fate of commutators 
(i.e. the quantum nature) in the scaling limit see the discussion in subsection 15.31 



5 Rigorous Results on the (Quantum) System 
in the Intermediate Regime and in the Scahng 



In this section we assume that the theory exists in the scaling limit provided 
that the scaling exponents have been appropriately chosen. Under this proviso 
we investigate its algebraic and dynamical limit structure. 

5.1 The Description of the System at Varying Scales 

In algebraic statistical mechanics we describe a system with the help of an observ- 
able algebra, A, a state, or expectation functional, < o >, a time evolution, 
at- Frequently one also employs the GiVS'-Hilbert space representation of the 
theory, introduced by Gelfand, Naimark, Segal (see e.g. [12]) • We already gave 
a brief discussion of these points in jS^. But as the approach of the scaling limit 
is quite subtle both physically and mathematically, we would like to give a more 
complete discussion of some of the topics in the following. 

We begin with fixing the notation and introducing some technical and concep- 
tual tools. Expectation values of elements of the underlying observable algebra, 
A, at scale "0", are given by 



where different indices may denote different elements of the algebra, different 
times etc. The dynamics is denoted by 



Limit 




(74) 



at{A) = A{t) or At , teR 



(75) 



space translations by 




(76) 
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at,M)^Mt,x) (77) 

Given such a structure, we can construct a corresponding Hilbert space rep- 
resentation (for convenience, we use the same symbols for the elements of the 
original algebra and their representations in the GiV5'-representation) . 

cu^n, uj{A{l) ■ --Ail)) = ifl\A{l) ■ ■■A{l)n)GNS (78) 

at^Ut, with at{A) ^UfA-U^t (79) 

The averaged or renormalized observables, A —>■ Ar, at scale R are a subset 
of elements contained in the original algebra, A. We denote the subalgebra, 
generated by these elements, by with Ar C A. We can decide to forget the 
finer algebra. A, and define the algebra on scale R: 

Definition 5.1 We define the system on scale R by 

u;(R)(A^R)):^u;{An) (80) 

af )(A(«)) := (a,(A))(«) (81) 

a^x\A^^'^) := (A(i?X))(^) (82) 

more specifically, we define the objects on the Ihs implicitly (via the GNS- 
reconstruction) by the following correspondence 

(A(«)(ti,Xi)---A(^)(t,,XO)(R) := {An{h,RX,)---An{ti,RXi)) (83) 

Remark: Note the different treatment of time and space-translations. We will 
come back to this point (which has remarkable physical consequences) below in 
connection with critical slowing down. 

Theorem 5.2 From the above we see that on each scale we have a new theory, 
S^^^ {S standing for "system"), which we get by reconstruction from the above 
hierarchy of correlation functions, in particular, a new, non-isomorphic algebra, 
A^^\ and a corresponding GA^yS'-Hilbert space representation. We emphasize 
that the coarse-grained dynamics is also physically different (despite the seeming 
similarity of the expressions on both sides of the above definitions) . 

If the scahng hmit does exist, we have, by the same token, a scahng limit 
system denoted by 

5°° = (a;°°,^~,Q;~,Q;^) (84) 

with 

{A'^{h,X,)---A^{ti,Xi)) = ^m{An{h,RX,)---AR{ti,RXi)) (85) 
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The proof is more or less obvious from what we have said above. 



□ 



Corollary 5.3 We generally assume that at is strongly continous on A. By the 
above identification process we can immediately infer that both al^"* and 
are also strongly continuous on the corresponding algebras, A^^\A°°. By the 
same token, we can infer that o;^^-' and u!°° are iiTMS'-states at the same inverse 
temperature /3. 

Remark 5.4 Strong continuity can be generally achieved by going over to smoothed 
observables, i.e., by averaging the observables with smooth functions of, say, com- 
pact support in the time variable. 

Proof of Corollary: Note that the original time evolution "commutes" with the 
scale transformation in the sense described above. This yields the mentioned 
result for all finite R. We have in particular that for suitable elements (for the 
technical details see [TBj) 

{B^^\t) ■ = (A(^) ■ i?(^)(t + zP))^n) (86) 

and there exists an analytic function, F^{z), in the strip {z = t + ir, < t < (3} 
with continuous boundary values at r = 0,/3: 

FiHit) = (A(^) ■ i?(^)(t))(K) , FS(t + ^P) = (A(^) ■ B(^\t + ^/3))(«) (87) 
This is equivalent to the following equation (cf. |16j): 

J o;(^)(A(^) ■ 5(^)(t)) • f{t)dt = J cu(^)(5(^)(t) ■ A(«)) ■ fit + if3)dt (88) 

for f E T). As f{t + ijS) is of strong decrease in t the limit R ^ oo can be 
performed under the integral and we get the same relation in the scaling limit. 
The above mentioned equivalence of this property with the ii'MS'-condition shows 
that the limit state is again KMS. This proves the statement. □ 

Remarks: i) Note what we have already said in 0. One reason for the non- 
equivalence of the algebras on different scales stems from the observation that, 
in general, 

Ar-Br^{A- B)r (89) 

Furthermore, in the scaling limit, many different observables of A converge to 
the same limit point, for example, all finite translates of a fixed observable, 
ii) A corresonding result in a slightly different context was also proved in 0]. 
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5.2 The Scaling Limit Theory as a Quantum Field Theory 

We have seen in sect. 2.3 that the scahng hmit of the correlation functions for 
the block spin observables is not fully scale invariant but only asymptotically so 
(while the short range details of the original microscopic correlations, encoded in 
the function F{xi — X2), have been integrated out, there remains an integrated 
effect of the initial block-function, f{x) ). 

This observation runs a little bit contrary to the general folklore, in which the 
various limit procedures are frequently interchanged and identified without full 
justification. We will exhibit the true connections between the various expressions 
in the following. 

With f{x) now being a general test function of e.g. compact support, we have 
from sect. 2.3, making now the dependence on / explicit 



R 

with 



hm {Arj{RX,) . Arj{RX2)) = canst ■ [ \y + y |-(-") • / * f{y)dy (90) 

Arj{RX) = i?-("+°)/2 . J A{RX + x) ■ f{x/R)dx (91) 
We now rewrite the limit correlation function as 

{Af{X^)-Af{X2))^ J {A^{x, + X,) ■ A'^{X2 + X2)) ■ f{xi)f{x2)dx,dx2 (92) 
that is, we identify 

Af{X) = J A°°{x + X)- f{x)dx (93) 

with A°°{x) now having rather the character of a field or operator valued distri- 
bution. 

We have that 

(i°°(Xi) • i°°(X2)) =: W°°{Xi - X2) = const ■ \Xi - Xsj-^""'^) (94) 

Corresponding results would hold for the higher correlation functions, that is, we 
arrive at 

Conclusion 5.5 In contrast to the block observables, AJ^, the field, A°°[x), 
displays the full scale invariance. 

The field, A'^{x), can, on the other hand, be directly constructed by means of 
a related limit procedure, which is however not of block variable type. We start 
instead with the unsmeared observables and take the scaling limit, R ^ 00 

lim{Aii{RXi) ■ Ar{RX2)) with Ar{RX) := • A{RX) (95) 
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and n — 7 = {n — a)/2. 

Remark: The extra scaling factor, i?", replaces the missing integration over the 
test function /ij, the support of which increases like i?". 

Performing the same calculations, we see that the above limit is equal to 
(A~(Xi) ■ A~(X2)). We hence have 

Conclusion 5.6 The fully scale invariant limit theory is achieved by taking the 
limits 

lim(i«(i?Xi) ■ ■ ■ Ar{RXi)) =■ W^{Xu . . . , (96) 

R 

The same construction holds of course for the intermediate scales; we define 
A'^^\X) by the following identification 

(i(«)(Xi) . . .i(«)(XO)(i?) := ■ {A{RX,) ■ --AiRXi)) (97) 

and have for the observables, A^ , defined above 

Af\x) = j A^^\X + x)f{x)dx (98) 
(which can e.g. be checked by direct calculation). 

5.3 The (Non)-Quantum Character in the ScaUng Limit 

In subsection B of section 3 of ;,8j, we already discussed the limiting behavior 
of commutators of scaled observables. In the regime of normal scaling, that is, 
scale dimension 7 = n/2, we found that commutators are non- vanishing in the 
generic case in the limit. This means that in general the resulting limit theory 
is non-ahelian (but quasi-freel) . Perhaps a little bit surprisingly, the situation 
changes at the critical point, where the scale-dimensions are, typically, greater 
than n/2 for at least some observables! 

We make the same observation as Sewell in namely, commutators of 
certain "critical" observables vanish in the scaling limit, i.e., the corresponding 
limit observables are loosing (at least) part of their quantum character . 

Remark: We think that the observation that fluctuations and critical behavior 
at the critical point are typically of a thermal and not of a quantum nature, does 
somehow belong to the general folklore in the field of critical phenomena, but 
we are not aware at the moment that this fact has been widely discussed in the 
literature in greater rigor. Some remarks can e.g. be found in connection with 
so-called (temperature-zero) quantum phase transitions in [T7] or [T3] and further 
references given there. 

On the other hand, related phenomena were observed in the context of spon- 
taneous symmetry breaking in sect. 6 of [H] and for certain models by Verbeure et 
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al in . A careful analysis of the behavior of commutators in a slightly different 
context can also be found in P3j . 

We have the following result. 

Theorem 5.7 Let A, B be strictly localized observables with 7a + 7b > n. We 
then have 

\x^\\[Ar,Br]\\=Q (99) 
Proof: With 7^ + 7^ > n we have 

\\[AR,Bn]\\<R-^^^^^^^ ■ j \\[A{xr),B{x2)]\\ - f{xr/R)f{x2/R)dx^dx2 

= ^-(7A+7B). j \\[A,B{y)]\\ - f{xjR)f{{x^ + y)/R)dx^dy (100) 

By assumption A, B have bounded supports, Va, Vb C so that 

[A, B{x)] = for + X n Va = (101) 

From the support assumption we immediately infer that the above double integral 
is actually a single integral as the commutator on the rhs vanishes outside a set, 
S, of finite diameter. We get 

hm II [Ar, Br] II < const' ■ R-^^^+^b) j ^^(^) . f^a:^/R)f{{x, + y/R))dx,dy 

= const' R-^^''^^''^ ■R'' j xsiy)-f{x'i)f{x'i+y/R)dx'^dy < const-lim/?"-(^-^+^^) = 

(102) 

as 7a + 7b > ■'^ by assumption. □ 

Corollary 5.8 We arrive at the same result if A, B are not strictly localized but 
fulfill a norm estimate of the form 

\\[A,B{y)]\\=:F{y)eL\R-) (103) 

Proof: We have 



^-(7a+7b) . I p^y^ . f^a:i/R)f{{xi + y)/R)dxidy = 

= R-^^^+^b) -R"- j F{p/R)f{p) ■ f{-p)dp (104) 
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We can again perform the i?-limit under the integral and get the hmit expression 



^n-(7A+7B) . ^(0) ■ j f{p) ■ fi-p)dp (105) 

for R —>■ oo. □ 
A simple example where different renormalisation exponents naturally arise is 
the following. Take a limit observable, A°°{X), and consider its spatial derivative, 
dxA°°{X). Then we have in a slightly sloppy notation (the limit being taken in 
the sense, described above): 

dxA'^{X) = \imdx{R~^^ ■ j A{x + RX) ■ //j(x)rf"x) 

= lim(i?(-^-^+^) ■ j {d^A){x + RX) ■ fR{x)rx) (106) 

That is, dxA = i[P,A] has to be scaled with a different scale exponent. Physi- 
cally, this can be understood as follows. With fni^x) = f{\x\/R) simulating the 
integration over a ball with radius R, a partial integration in the above formula 
shifts the dr^ to the test function, fni^x). As d^fni^x) = in the interior of the 
ball, the averaging goes roughly only over the sphere of radius R instead of the 
full ball. This has to be compensated by a weaker renormalisation. 

Another result in this direction can be found in [Sj sect. 6, in connection with 
the canonical Goldstone pair in the context of spontaneous symmetry breaking. 

Further possible candidates are the time derivatives of observables as, for 
example, in (AA). Fourier transformation yields an additional prefactor, uj"^ in 
the spectral weight, Faa{^^, k). The ii'MS'-condition leads to another constraint: 

Fab{uj, k) = {l- e-^'^)-i ■ CiA,B]i^, k) (107) 

A combination of such properties shows, that in the scaling limit, the vicinity of 
(a;, /c) = (0,0) is important (see the discussion in section Hj). 

From covariance properties (as e.g. in models of relativistic quantum field 
theory) one can infer additional information about certain characteristics of the 
energy-momentum spectrum. For arbitrary models of non-relativistic many-body 
theory, however, the situation is less generic and typically model dependent. 

Remark: We had several discussions with D.Buchholz about this point, which 
are gratefully acknowledged. This applies also to the following subsection. 



5.4 The Nature of the Limit Time Evolution and the 
Phenomenon of Critical Slowing Down 

We argued above that the appropriate choice of the respective scaling dimensions 
of the observables under discussion is a subtle point and perhaps, to some extent. 
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even a matter of convenience. After all, one may have some freedom in the choice 
of the subset of observables which is to survive the renormalisation process. 

We will not give a complete analysis of all possibilities in the following but 
rather emphasize one, as we think, particularly remarkable phenomenon, namely, 
the phenomenon of critical slowing down. As in the preceding discussion, we 
choose two observables. A, B, with 7a + 7_b > implying that the limit commu- 
tator vanishes. We assume this also to hold for non-equal times (which follows 
from a Li-cluster condition as in the preceding subsection), at least on the level 
of two-point functions, i.e. 

([A-,S-(t)])oo = (108) 

and get the following theorem: 

Theorem 5.9 Under the assumptions being made, we have 

{A°° ■ B°°{t))^ = const for aWt eR (109) 

Proof: As the limit state is again a i^'MS'-state, the vanishing of the above 
commutator implies that the analytic function, F^^{z), fulfills 

FTsit) = F^^it + (110) 

for all t. F^q{z) can hence be analytically continued to the whole C-plane and is, 
furthermore, a globally bounded analytic function, hence a constant by standard 
reasoning. □ 

We see that the subclass of limit observables, which has vanishing limit com- 
mutators (see the preceding subsection), has, by the same token, time indepen- 
dent limit correlation functions. As these pair-correlation functions are usually 
connected with characteristic observable properties of the system (generalized 
suszebtibilities, transport coefficients etc.), this has remarkable physical conse- 
quences. The corresponding phenomenon is called critical slowing down. For 
a review of the physical phenomena see e.g. In physical terms, the phe- 

nomenon can be understood as follows. 

In the critical regime, the patches of strongly correlated degrees of freedom 
become very large and extend practically over all scales. That is, a reorientation 
of such clusters or a response to external perturbations takes, if viewed on the 
microscopic time scale, a very long time. In the scaling limit this time scale goes 
to infinity. If one wants to see observable dynamical effects on the macroscopic 
level one must hence scale the time variable also. For the unsealed time we have 
in the limit i? — *• cxo: 

d/dt{A°° ■ B°°{t))^ = (111) 

This is the same as 

■ [H^,B^{t)])^ = \im{AR ■ [H,BR{t)]) = \imd/dt{AR ■ BR{t)) (112) 
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(cf. subsection 15.11 H is the microscopic Hamiltonian) . 

What one now has to do is obvious. We have to compensate the vanishing 
of the above expression in the hmit by inserting an appropriate scale factor in 
the time coordinate. Instead of B{t) we take B{R^ ■ t) with 5 so chosen that 
the hmit expression is non-vanishing. Note that differentiation with respect to 
t now yields an explicit prefactor R^. This fixes the macroscopic time scale, tm, 
for these processes. We define 

■ B^{tJ)^ := lim{An ■ Br{R' ■ tj) (113) 

-ft 

Physically the effect can be understood by inspecting the middle part of equation 
flll2p . The support of Bji(t) spreads with time. This spread is more pronounced 
if we take R^ ■ t instead of t. By the same token the overlap with the Hamiltonian 
(which is basically translation invariant) increases with R ^ oo while t is kept 
fixed, thus yielding the non- vanishing limit. 

It may happen that other observables may evolve on different macroscopic 
time scales so that the construction of a coherent common macroscopic limit 
time evolution may not be straightforward. Such more detailed questions have 
to be separately studied for the various model classes. As we have studied a 
concrete model class in section |3J we can make much more precise statements if 
we have some information about the energy-momentum spectrum in the vicinity 
of (cj, k) = (0, 0). In that section we arrived at the following results: 

F^^iu, k) = il- exp{-f3Lu))-' ■ c^b(^, k) (114) 

and 

F,i^g{uj, k) := 4{k) ■ 6{uj - {e{k) - /x)) + Jl{k) ■ 6{uj + {e{k) - /i)) (115) 

with 



4{k) = (1 - exp i-mk) - ■ m (116) 



and 

For e.g. 



Jt{k) = exp {-P{e{k) - ii)) ■ Jl{k) (117) 



{AR{RXi,t)-AniRX2,0)) (118) 
we have to study expressions like 

^-2(7^-n/2). f e^(Kfc/«)-M)*.(i_exp(-/3(e(A;/i?)-^)))-i-j(fc/i?)-|/pd"A; (119) 



Again, the situation is normal for /i 7^ but becomes singular for /i = (cf. 
section HI) . 
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We concentrate on the case e{k) ~ \k\°' for \k\ — > and j{k) continuous 
and 7^ in A; = (cf. sectiorH] assumption |3] 6). In order to have a non- vanishing 
hmit correlation function we have to choose 

7A = {n + a)/2 (120) 

In the case where j{k) vanishes polynomially in A; = we have to make a corre- 
sponding choice, as has been described in the mentioned section. 

If, furthermore, we want to have a non-trivial time evolution in the limit 
R ^ (yo, we have to scale the microscopic time like 

t = R'-T so that |A;/i?|" -t = ■ r (121) 



Conclusion 5.10 If we are in the situation, described in section |3J having a 
singular spectral contributution with quasi-particle-like dispersion law e{k) ~ 
near = 0, we have to scale the microscopic time, t, like t = R" ■ r, in order to 
arrive at a non-trivial limit time evolution in the variable r. 

Acknowledgement: Several discussions with D.Buchholz are greatefully ac- 
knowledged (see also the remark at the end of subsection 5.3). 
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